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the derivatives of the system output go to zero at the instant that the 

output variahle arrives at a steady-state value. 

The purpose of this study is to apply the method presented in ref- 
erence 7 for determining the optimum switching criterion to a third- 
order system having complex conjugate roots. Furthermore, this paper 
presents a relatively simple analog-computer method for obtaining this 
criterion and the restilting optimum transient responses. 

In order to introduce the methods involved, the switching criterion 
is first determined for a second-order system that is descriptive of an 
ideal missile roll control system. The methods are then applied to a 
third-order system that is descriptive of an ideal missile acceleration 
control system. The application of this method to this system is com- 
plicated by the presence of complex conjugate roots. 

Transient responses of a third-order limited-linear system are also 
included in this paper. These transient responses present a basis of 
comparison for the optimum third-order contactor control systems. 


SYMBOLS 


constants of integration 


b 

C 

c 


E 

Ei 


constant coefficient in airframe normal -acceleration 
transfer function, radians/sec 

nondimensionalized constant, c/b^ 

constant coefficient in airframe normal -acceleration 
transfer function, (radians/sec) ^ 

nondimensionalized system error 

nondimensionalized system error input 


g acceleration due to gravity 

j = 

velocity constant of servomotor in normal -acceleration 
control system, deg/sec/volt 

K 2 proportionality constant in limited-linear normal - 

acceleration control system, volts /g unit 
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K3 


K4 


K5 


n 




S 

T 

t 

^na 




Ar 

A = 
6 


e 


Ai, A2^ 


T 


feedback gain constant In limited-linear normal- 
acceleration control system, volts /g unit/sec^ 

velocity constant of airframe transfer function in 
roll control system, deg/sec/deg 

static gain constant of airframe transfer function 
in noimial -acceleration control system, g units/deg 

control-system normal -acceleration output signal, 
g units 

control-system normal-acceleration input signal, 
g units 

Laplace transform variable 

roll-control-system or airframe time constant, sec 
time, sec 

coordinate of w phase space 
locus of switching points in w^ space 
forcing function for mathematical analysis 
actuating signal for servomotor, volts 


control-surface deflection, deg 

control-system actuating signal, g units or deg 

characteristic roots of normal-acceleration control 
system 

dimensionless time 

control-system roll-angle output signal, deg 




control-system roll-angle input signal, deg 
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A dot over a quiantity denotes the derivative with respect to time t. 

A primed quantity denotes the derivative with respect to dimensionless 
time T. 

The subscript a is 
The subscript n is 

DEFINITIONS 


employed as a switching index (a = 0, 1, 2, and 5) • 
eii 5 )loyed as a coordinate index (n = 1, 2., and 5) • 


Contactor: A nonlinear element which has an output constant in magni- 

tude with the sign of the input and, furthermore, has zero output when 
the input is zero. 

Optimum transient response: A transient response that reaches the steady- 

state value in a minimum of time with no overshoot. 

Switching criterion: The system conditions at which contactor reversal 

must take place to yield optimum transient responses. 

Locus of first switching points: The locus of system coordinates required 

to define the first contactor reversal points in an optimum contactor 
control system. 

Locus of second switching points: The locus of system coordinates 

required to define the second contactor reversal points in an optimum 
third-order contactor control system. 

Zero trajectories; The phase-space trajectories passing through the 
origin. 


SYSTEM DESCRIPTION 
Second-Order System 


A block diagram of the first system considered in this investiga- 
tion is given in figure l(a). This system is descriptive of a missile 
contactor roll control system. A step input 0^ to the system causes 
an error e to be applied to the contactor element. The sign of this 
error signal actuates the contactor which then applies a step aileron 
deflection 6 to the control surface. A reversal of the sign of the 
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error signal when the error reaches zero causes a corresponding reversal 
of the contactor output. In order to facilitate the analysis, 6 is 
limited to ±1°. The magnitude of the aileron deflection is Included in 
the airframe velocity constant. The performance of this system is char- 
acterized by a large initial overshoot and a decaying oscillation. 

The contactor system of figure l(a) was then modified by the inclu- 
sion of a computer element controlling the contactor input to cause an 
optimum system respnse. Figure l(b) is a block diagram of the modified 
system. In response to a step input this element actuates the con- 

tactor which applies a step aileron deflection 6 to the airframe. 

Unlike the system of figure 1(a), the sign of 6 is reversed before the 
error reaches zero. In the optimum case, the sign reversal takes place 
at the appropriate time to cause the error and the first time derivative 
of the error to reach zero simultaneously. The return of the control- 
surface servo to a neutral position removes the control- surface deflec- 
tion, and the system remains at rest. The neutral relay position required 
to keep the system at rest is not shown in the block diagrams j therefore, 
figures l(b) and 2 apply to the systems only when they are not at this 
rest or steady-state position. 


Third-Order System 

A block diagram of the optimum third-order contactor control system 
considered for this study is given in figure 2. The performance of this 
system is similar to that of the optimum second-order system previously 
described with the exception that, because of the order of the system, 
the coxiiputer element must actuate the contactor in such a manner that 
the error and its first two time derivatives go to zero simultaneously. 
In order to accomplish this, a minimum of two switchings is reqiiired. 

In response to a step input of acceleration, the contactor element of 
the optimum third-order system applies a constant rate of pitch control- 
surface deflection to the airframe. In a manner similar to the second - 
order system, the contactor output is restricted to il by including the 
actual magnitude of this output in the constant K^. The resulting rate 

of control -sxu-face deflection 6 causes the error to decrease at an 
increasing rate. Then, the contactor output is reversed. The error 
continues to decrease but with a decreasing rate. Once more reversal 
of the contactor output causes the error, the error rate, and error 
acceleration to go to zero simultaneously. Figure 3 illustrates a typi- 
cal response of this system to a step input of the controlled variable. 
This figure contains time histories of e, k, e, A^., and n^. 
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Figure 4 is a block diagram of a limited -linear third-order normal- 
acceleration control system. The forward loop gain K£ a-nd the feedback 

pt 

gain Kj of this system were adjusted to minimize 
5g step input. 


^0 


|e (t) I dt for a 


METHOD OF AMLYSIS 


The analysis is concerned with the determination of the switching 
criterion necessary for the optimum performance of a contactor control 
system. The mathematical analysis is based primarily on that developed 
by Bogner in reference 7. The transfer-function coefficients used for 
the systems discussed are presented in table I. 


Second-Order System 

The equation for the simple contactor roll control system shown in 
figure l(a) can be written as 


T ^ M = K4& (1) 

dt^ 

Inasmuch as € = if only step inputs of position are considered, 

equation (l) can be written: 


T 




( 2 ) 


Substituting e 



and performing the integration yields 


i e + € - K 4 S loge(e + K 45 ) = Bq 


(3) 


where Bq is the constant of integration. Equation (j) yields the phase 

portrait for the system from which any of the possible system trajectories 
may be obtained. For the system of figure l(a), the sign of 6 reverses 


MCA TN 57^4-5 


7 


as e changes sign. A typical phase-plane plot of the response of the 
system to a step input of position obtained from an analog computer is 
shown in figure 5^ and a transient response for this system is presented 
in figure 6. 

Evaliiation of the constant of integration Bq in equation (5) for 
the phase-plane trajectories through the origin gives 

i 6 + I - Ki,6 + 1^ = 0 (4) 

Because of the dual value of 6 (that is, 6 = tl°), equation (4) 
is actually two equations: 

I e + 6 - Ki, log^(l + i^) = 0 (5a) 


i e + I + Ki^ logg ^1 - i-j = 0 (5b) 

Inasmuch as there exists only one solution of the equation ^ = f(e,e) 

through each point in the phase plane, equations (5a) and (5b) represent 
the only paths along which the representative point of the system may 
enter or leave the origin. These are called the zero trajectories. 

Figure 7 shows zero trajectories obtained from equations (5a) and 
(5b) on an analog computer and by direct nimierical substitution. Por- 
tions of these zero trajectories will later be used for the optimum sys- 
tem of figure l(b). If equation (2) is set up on an analog computer, a 

phase-plane plot of e against e obtained as the representative point 
moves from the origin with 6 = +1° will give curve oa of figure 7. Simi- 
larly, a solution of the same equation with 6 = -1® yields curve ob. 

The substitution of -t for t in eqmtion (2) results in the following 

equation: 
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A phase -plane plot of e against e obtained for equation (6) in a 
similar manner from an analog computer with 5 = +1° and 5 = -1° yields 
the curves oc and od, respectively, of figure 7. (See ref. 8.) 

The segments of the zero trajectory od and oc entering the origin 
can be utilized as the contactor switching cixrve. If the sign of 6 is 
reversed at the instant a system trajectory first intersects a segment 
of the zero trajectory entering the origin (at a point such as (a) in 
fig. 5), the system trajectory will then follow the zero trajectory to 
the origin and execute a deadbeat response. This switching criterion 
will yield an optimum transient response. 

The block diagram of the second-order system, modified by the addi- 
tion of a computer element to reverse the sign of 5 when the represent- 
ative point first intersects the switching curve, is presented in fig- 
iire l(b). This computer element senses the region of flg\rre 7 in which 
the representative point of the system trajectory lies, and produces the 
proper contactor output. Whenever the representative point crosses the 
segment of the zero trajectory leading to the origin, the contactor out- 
put reverses. Finally, at the instant this point enters the origin, the 
contactor output goes to zero. Figure 7 also illustrates typical phase- 
plane responses of this system to step inputs of position, and figure 8, 
a number of transient responses obtained on an analog computer. 


Third -Order System 

In the second-order contactor system, the switching criterion causes 
e and e to go to zero simultaneously. For an optimum third-order sys- 
tem, it can be shown that £, €, and € must go to zero simultaneously. 

It is shown in reference 7 that two switchings are required to do this. 

In a third-order system the second or final switching c\rrve in phase space 
leads to the origin and, hence, is still called the zero trajectory. In 
general, the locus of first switching points in a third-order system con- 
sists of a surface containing the zero trajectory. If, however, the 
Inputs to the system are restricted to steps of the controlled variable, 
the locus of first switching points becomes a line lying in this surface. 

In the following analysis the equation for the third-order system 
is presented. With the inputs restricted to steps of the controlled 
variable, a pair of simultaneous equations for the locus of first switching 
points is obtained. The zero trajectory is also developed from the equa- 
tion for the system. An analog-computer method for obtaining this 
switching criterion is introduced and used to check the analytically 
derived criterion. Transient responses are then obtained from the con- 
tactor system using this switching criterion. 
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Analytical method .- A block diagram of the third -order contactor 
acceleration control system studied is given in figure 2. The differen- 
tial equation of the basic system in terms of the error variable can be 
written as follows: 


•e(t) + be(t) + c€(t) = K1K5C 


(Ta) 


In the appendix, where a complete derivation of the switching criterion 
is presented, equation (7a) is nondimensionalized and written as 


E’" (t) + E"(t) + CE'(t) = A 


(Tb) 


Equation (7b) is in general form inasmuch as time and amplitude scales 
may be adjusted for any desired \mdeunped natural frequency and system 
velocity constant. The resialts obtained can be applied to any similar 
system with the same damping ratio (0.l4). 


In the following analysis, after two transformations of variables, 
the switching criterion is determined in terms of Wj_, W2, and w^. 

These transformations were used to facilitate the analysis. The varia- 
bles W]_, W2^ and wj are defined by the following three eqviations: 


w-i (t) = E(t) + — — E' (t) + — i — E" (t) 


^ 2^5 


A2A5 


W2(t) = - 


A2A5 


E'(t) 


A2A5 


E"(t) 


(8) 




Where the characteristic roots of the system are represented by A. 

In determining the loci of switching points, the variables Wj^g^ 

are introduced. The locus of first switching points will be written 
in terms of wn, W2p^ and wji and the zero trajectory, or locus of 
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second switching points, in terms of '^22> ^J>2 hoth 

loci lying in the W]^, W 2 , and w^ phase space. It is shown in the 

appendix that the locus of first switching points can be represented by 
the following two equations: 


1 - J»5i) . 


1 - 


^2^ (^3 - >>2) 


2 An 


(^21 + J^3l) 


A3/A2 


(9) 


and 


“1 


1 - 


^2^ (^3 “ ^2) f ^2^^3l 

i-2 exp Wn n ^ 

2 Ai V 


W21 + JW52. + 


2 An 


"^2^ (^3 - ^2) 


1 ± 




(^3 ~ ^2) 


2 Ai 


exp wii 




2^3 

Ai 


w, 


21 - JW5I 


2 An 




^3^ (^3 ■ ^2) 


A2/A3 


( 10 ) 


Figlire 9 contains a plot in the W5W2 plane of values of W21 and 
w^l which satisfy equation (9)* Inasmuch as equations (9) and (lO) are 
simultaneous equations^ the values of W21 and w^^ which satisfy equa- 
tion (9) were substituted point by point into equation (lO), and the 
resulting equation solved for w^^q^. The values of “thus obtained 

are shown plotted against ^21 '^51 figures 10 and 11. This 

numerical evaluation was done with a digital computer. Figure 12 is a 
sketch of the locus of first switching points in the upper half of the 
w phase space . 


It is further shown in the appendix that the locus of second 
switching points or zero trajectory can be expressed by the following 
pair of simultaneous equations: 


^2^ + 


2) _ j ^2^^3^(a3 - A2) ^ 
/ Aq 


'32 = ^3^ 


exp 



w 


12 



i'^2^ 

Al 



exp 


(11) 
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“ A2) 


W22 



exp 


^ Al 




exp 




(12) 


The W]_Wj projection of the locus of second switching points is plotted 

in figure 13* Because the inputs are restricted to steps of the con- 
trolled variable^ only the initial portion of this trajectory on either 
side of the origin is used. 

Analog-conrputer method . - The two loci of switching points can be 
obtained by an analog-computer method. The determination of the switching 
criterion by this method involves the use of the original equation for 
the system (eq. (7t>)) on an analog coir^uter. Inasmuch as this equation 
is used, the criterion can be obtained in terms of the original system 
variables. However, for the presentation given here, the transformations 
of equations (8) were incorporated with equation (7b) in the analog- 
computer setup. This permitted the use of the analog-computer method 
to check the analytical results. 


Transformed system coordinates: The actual determination of the 

switching criterion by an analog-computer method is accomplished in 
three steps. First, the zero trajectory or locus of second switching 
points is obtained. Second, the w^W2 projection of the locus of the 
first switching points is determined. Finally, after the complete space 
trajectory of the locus of second switching points and the w^W2 pro- 
jection of the locus of first switching points are detennined, the 
remaining W]_ coordinates of the locus of first switching points are 
obtained. In detail, this procedure is as follows: 

(1) A solution is obtained in negative time from the origin of the 
w phase space with *i:A. The resulting space curve is the path along 
which all trajectories entering the origin must travel. As in the second- 
order system, this space curve is called the zero trajectory. 

(2) For step inputs (w2_,0,o) to the system in positive time, all 
initial trajectories have the same w^W2 projection. Any subsequent 

switching from such initial trajectories must also have this projection. 
Hence, the w^W2 projection of the locus of first switching points is 

determined by operating the system in positive time with tA from the 
origin or any point on the w^ axis and recording the w^W2 projection 

of the resulting trajectory. 
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(5) The zero trajectory and the wjW2 projection of the locus of 
first switching points have been determined. The W]_ values of the 
locus of second switching points remain to be foiind. Figure l4 shows 
the vy/f2 projection of these two loci. If the system leaves the ori- 
gin (point o of fig. l4) in negative time along the zero trajectory and 
A is switched at some point b between the origin and point a of curve oa, 
the value of w^ recorded at the intersection of the resulting trajectory 

6ind the w^W2 projection of the locus of first switching points (points c, 

d, e, and f of fig. l4) is the remaining coordinate of this space curve. 

If the sign of A is again reversed at this Intersection, the trajectory 
will follow the locus of first switching points to the origin. By varying 
the point on the zero trajectory where the system is first switched in 
negative time, the space curve of the locus of first switching points can 
be defined. 

Values of the zero trajectory obtained by the method discussed are 
plotted in figure 15 along with solutions of the analytical expression. 
Projections of the locus of first switching points obtained with this 
method are shown with digital -computer solutions of the analytical 
expressions in figures 9, 10, and 11. It was found, during the analog- 
computer determination of the locus of first switching points, that the 
W]^ values associated with this space curve must be determined with care. 
Because of the switchings Involved, it is possible to introduce consider- 
able error at this point. 

The phase-space operation of the system in response to a step input 
of the controlled variable is shown in figures 15 and 16. The phase - 
space trajectory begins with the initial value of error at point a. It 
travels along curve ab of these two figures until it intersects the locus 
of first switching points at b. At this point, the sign of Aj. is 

reversed. This reversal changes the direction of the trajectory, which 
now travels along curve be until the zero trajectory is reached at 
point c. Here the sign of is again reversed and the system travels 

along the zero trajectory to the origin at point d. 

Original-system coordinates : The analog method can be s-ummarlzed 

in terms of its use to determine the switching criterion in terms of an 
original system error variable. The equation for the system is set up 
on the computer and a solution for £(t) is obtained in negative time. 

The trajectory e(-r) = f [i(-T), e(-T), obtained in this manner is 

the system locus of second switching points or zero trajectory. The 
equation of the system is then set up in positive time and the projec- 
tion of the trajectory in the e(t)e(t) plane is recorded. This pro- 
jection coincides with that of the locus of the first switching points 
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in this plane. The remaining coordinate e(t) is obtained by plotting 
the system trajectory on the e(t)e(t) plane as the system leaves the 
origin in negative time and then reversing the sign of A until the 
trajectory intercepts the previously mentioned projection plotted in pos- 
itive time. The value of €(t) at this Intercept is the remaining 
required coordinate of the locus of first switching points . The proce- 
dure is repeated until sufficient points on this locus are established. 


DISCUSSION OF THIED-ORDER TRANSIENT RESPONSES 


Transient responses obtained from the optimum contactor normal- 
acceleration control system for several step -input amplitudes are shown 
in figure IJ. These responses are evaluated by comparison with those of 
an equivalent limited-linear system. The linear system utilized for this 
purpose is shown in the block diagrajn of figure 2 . The maximum rate of 
control-surface deflection of the linear system was limited to that of the 
relay system. In addition, the forward loop gain K2 and the feedback 


for a step input of 5 g» A series of transient responses obtained from 
this system for several step-input amplitudes are shown in figure 18. 

A comparison of figures I7 and 18 shows that in the optimum system over- 
shoot is eliminated. In addition, this comparison also shows that there 
is a considerable superiority in rise time of the optimum system for 
input amplitudes below those at which the linear-system gains are adjusted. 
For larger input amplitudes, this superiority is less because the limited- 
linear system operates at its velocity saturation limit for a large per- 
centage of the transient time. For the optimum case, the control surface, 
in response to a positive unit step input n^, initially travels at a pos- 
itive constant rate and then at the appropriate time travels at a negative 
rate. At the second switching point, the control surface travels at the 
positive constant rate again until the output n arrives at the steady- 
state value. This control-surface response is typical of the optimum 
contactor control system and is proportional to the integral of Aj.(t) 
of figure 3» 


This paper presents a switching criterion that yields an optimum 
transient response to step inputs of the controlled variable for a third- 
order contactor system with complex roots. This optimum response has no 
overshoot and arrives at a steady-state value in a minimum of time in 
response to a step input of the controlled variable. 


gain Kj of the linear system were adjusted to minimize 



0 


CONCLUDING REMARKS 


14 


MCA TN 5745 


Two methods of obtaining this switching criterion, a direct mathe- 
matical method and the analog-computer method, are presented. In the 
direct mathematical method, two transformations of variables were applied 
to the original nondimens ionali zed equation of the system to determine 
the equations of the two loci of switching points in the phase space. 

The equations of these loci are in general form inasmuch as the time 
scale of the transient outputs as well as the amplitude scale may be 
adjusted for any desired undamped natxaral frequency and system velocity 
constant for the given value of dancing ratio. A digital computer was 
utilized to obtain the actual coordinates of the switching loci after 
the equations of these loci were established. 

The analog-computer method, although lacking the accuracy of the 
mathematical method, offers a relatively simple means of obtaining the 
two loci of switching points . In addition, these loci can be determined 
in terms of the original system error variable with no transformations . 

The transient responses of the optimum system were compared with 
those of the equivalent limited-linear system. The optimum contactor 
acceleration control system is superior with respect to the time required 
to reach steady-state values; however, for large step inputs of the con- 
trolled variable, this time advantage is less because the limited-linear 
system operates at the velocity saturation limit for a large percentage 
of the transient time. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 10, 1956. 
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APPENDIX 

DERIVATION OF THIRD-ORDER SWITCHING CRITERION 

Symbols 


G] 






ma-brix 


matrix 


0 10 
0 0 1 
0 -C -1 

o' 

0 
A 


transformation matrix 


inverse matrix P 


coordinates in space 




dr 


n 


matrix 


U2 

^5 


(n = 1, 2, and 5) 


'n 


principal coordinates of system (n = 1, 2, and 5) 


matrix 


a 

V2 

^5 




i6 
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Third-Order Switching Criterion 

The basic form of the equation for the third-order system of fig- 
ure 2 is 


JL = (Al) 

^ S(S2 + bS + c) 

When e = ni - n and inputs are restricted to steps of n^ so 
that Hi = iii = n’i = equation (Al) can be written as follows: 

e'(t) + be(t) + ce(t) = -KiK^c (A2) 

This equation can be simplified with the following substitutions : 

A = -Aj, 


T = bt 


E(t) = — — e(t) 
K^K^c 


C 


c 



This simplification results in 

E"'(t) + E"(t) + CE'(t) = A (A5) 


Let 

U]L = E(t) 
U2 = E' (t) 
U5 = E"(t) 


(A4) 


3L 
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From equations (A5) and (a 4) the following equations can be 


du-j 

dr 


= 


■2 _ 


= U7 


dU; 
dr 

dux 

= -CUp - Ut; + A 

dT ^ ^ 


In matrix notation, equations (A5) become 


dr 


^1 


0 

1 

0 


^1 


- - 
0 

U2 

= 

0 

0 

1 


U2 

+ 

0 



0 

-c 

-1 


"5 


A 


or 


i u] . A u] + g] 


The roots of the characteristic equation are 


^1 
A2 = 

A:z; = 


= 0 


-1 + 


^1 - 4c 


-1 - ^1 - 4C 


When 


jp = ^1 - 4C^ these roots can be written as follows; 


^1 
A2 = 

Ax = 


= 0 


-1 + jP 
2 

-1 - jP 


IT 

written : 
(A5) 


(a6) 

(at) 


(A8) 


2 
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Introducing a new coordinate space with the transformation 

matrix P so that 

u] = P v] (A9) 

permits writing the system equation (AT) in terms of the new coordinates 
as follows : 


A v] = P"^AP v] + P“1 g] (aid) 

dr 


When the principal-coordinate transformation is applied, the 
matrix P is determined so that the resulting P-lAP matrix is of the 
diagonal form. This selection results in (see ref. 7) 


P = 


111 
0 ^2 

0 >v2^ ^5^ 


From equation (All), where - ^^2 + = Ij 


D -1 - 


A5A2 


0 


A 5 A 2 

-1 


1 


^(^5 - A2 ) A 2 ( Aj - A2 ) 

-^2 1 

A3 (A3 - A2) A3 (A3 - A2) 


(All) 


(A12) 


After substitution of the values of eqmtions (All) and (A12) for 
the transformation matrix and its inverse, the matrix P'^AP can be 
evaluated as 
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P-^KP = 


0 0 0 
0 A2 0 

0 0 Aj 


(A13) 


The substitution of equations (A12) and (A13) into equation (AlO) gives 



^1 


000 




1 

1 1 
A^A 2 ^ 5^2 


"0 




0 

ro 

0 


V 2 

4. 

0 

I 

H 


0 

dr 

^2 



T 


^ 2(^5 - ^2) A 2 ^A 5 - 

A 2 ) 


t/"— 


0 0 A 5 


“XT-. 


0 

-A 2 1 




"^5 

l_ _ 



^5_ 


r< 

1 




which is equivalent to 


^ A 

d.T ^5^2 


(Al4) 


dv2 

dr 


A 2 V 2 


^2(^5 - A2) 


(A 15 ) 


dvj 

dt 


A 5 V 5 + 


A 

^5(^5 " ^2) 


(AI 6 ) 


Eliminating dr between equations (A14) and (AI 5 ) s-iid between equa- 
tions (Alij-) and (AI 6 ) and integrating results in the following eqmtlons : 


vi = 


A2^A5 


loge 


A 2 V 2 - 


^2(^5 - ^2) 


+ Bi 


(AI 7 ) 
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Vi = 




log<= 


Xjv, + 


A5(X5 - X 2 ) 


+ B 2 


(A18) 


Equations (AIT) and (AI 8 ) can be solved to yield a zero trajectory 
in terms of v^^, V 2 ^ and v^ and a locus of first switching points in 

terms of these same variables. These variables are, however, complex. 
In order to reduce these equations to real variables, it is first noted 
that equation (A9) when solved for v^ yields 


v] = p-1 u] (AI 9 ) 

The substitution of equation (A12) into equation (AI 9 ) leads to the 
following set of equations: 


V^^ = Uq_ + 


— i — Up + — — Un. 
A 2 A 5 ^ 2^5 ^ 


Vo = 


A 5 


A2(A3 - A 2 ) ^ A2(A5 - A 2 ) 


U:! 


A 2 1 

= : ^ — r 1^2 + •; — ^ 


^5(^5 - A2) ^3(^3 " ^2) 


(A20) 


The substitution of equations (A8) into equations (A20) gives 


^1 = ^1 + ^ ^2 + 


U7 


p ^ p 9 

+ 1 P^ + 1 


V2 = - — U2 + a 


. (1 - 3^) 


U 2 - Uj + j 


U^ 


32 + 1 3(32 + 1) 32 + 1 ^ 3(32 + 1) 


''5 - - “2 - J ^ U, - , uj 


(A21) 
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Inasmuch as V 2 and are con^jlex conjugates and V]^ is a real 

variable, the following change of variables will yield equations in terms 
of real variables : 


wi = vi 

W2 = V2 + vj 

^2 - V 3 


wj = 


J 


(A22) 


Substituting equations (A22) into equations (A21) yields 


4 4 

^1 = ^1 + — ^2 + — U5 

+ 1 + 1 


W 2 = 


4 4 

U 2 - U 5 


p2 + 1 " p2 + 1 


w, = ~ U 2 + 

p(p 2 + 1 ) p(p 2 + 1 ) 


(A25) 


or, in terms of A, 


Wj^ = Uj^ + — — U 2 + — U7 


^2^5 ^2^5 


Wp — — — — up — U3 

A2A5 A2A3 ^ 


^5^ + A 2 ^ 

w^ = -J ^ Up - j 


A2A3(Aj - A 2 ) ^2^3 (^3 “ A 2 ) 


(A24) 
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Substituting equations (A22) into equations (AIT) and (A18) yields 


wi = 




loSe 


^ '”' 5 ) ■ X2(^5^ X 2 ) 


+ Bi 


(A25) 


Wn = 


A 2 A 5 


5 logs 


'H^2 - J^5) + 




1(^3 - A2) 


+ B 2 


These two equations can be rewritten as follows : 

2A 


W2 + Jwj 2 


A2 ^(Aj - A 2 ) 


explW]^ 





^2^^3^ = 


B7 


(A26) 


(A2T) 


W2 


2A 

'"5 * . X 2 ) _ 


explwp 


^ 2 ^ 3 ^ 


^exp(-B2^ 


= B4 


(A28) 


A typical response of the system considered to a step input ^ of the 
controlled variable will begin at E(t) = E*(t) = 0, and E (t) = 0 

and terminate at E(t) = 0^ E*(t) = 0^ and E (t) = 0 in the E(t) 
phase plane. The initial point of this response will be called 0, the 
first switching point 1, the second switching point 2, and the terminal 
point or origin 5 . By using the transformations of equations (A24) and 


w, 


na^ 


the notation that Wj^ at a point a along a trajectory becomes 
this response can be shown to begin at wp = Ej_ = wpo> ^2 = 0, w^ = 0 
and terminate at wi = 0, W2 = 0, and Wj = 0 in the w phase space. 
This response is shown schematically in the following diagram. 


NACA TN 57^5 


25 



The first trajectory of this diagram is described by eqmtions (A2'j) 
and (A28) with substituted for A. The constants and have 

the same respective values at any point on this trajectory. If Bj is 
evaluated from eqmtion (A2j) at point 0 and equated to Bj eval\iated at 
point 1 on this trajectory, the following equation can be written: 


w. 


21 


2Ai 

* '“5^ ■ ) 


exp w^i 


^2 ^5 


Ai 


2Ai 

^22 (X5 - Xs) 




The application of this procedure to the three trajectories for both 
and results in the following six equations in terms of the ini- 

tial and terminal points and the two switching points : 
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From initial trajectory. 


W21 + jw^i - 


2Ai 


2Ai 


^2^ (^5 ” ^2 ^ ^2^ (Aj "42^ 


exp wii 






exp w^o 


^2^ 


2^ 


2Ai 


^21 - JW31 + _ 7,2) ^52(^5 - A2) 


exp Wii 


^2^3^ 

Ai 


exp Iwio 


^2^3^ 


from second trajectory, 


2A-, 


-"x,j 

/ A2^A7 

Ki 


2Ai 


“22 * ^"52 - ^ 2 ) 

A2^Aj\ 


exp 1 -w. 


12 




W21 - jvjx 


2Ai 


2Ai 


Aj2 - A 2 ) 


W 22 " j''^32 " 


^3^ (^5 ■ ^2) 


AoA- 


exp 1 


2^3 

Ai 


AoA 


exp l-w- 


12 


2^3 


and from zero trajectory, 

■W22 + J^32 

exp Wx2 


2An 


A2^(A3 - A 2 ) ^ _ 2Ai 


^^3 

Al 


^2^ (^3 " ^2) 


(A29) 


(A3O) 


(A3I) 


(A32) 


(A35) 
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W22 “ J ^32 


exp 



2^1 

^3 ^ (^3 ~ ^2) 

>^2^3^ \ 


2Ai 

^3^ (^5 ■ ^2) 


(A3i^) 


The elimination of ^12 •> ^22-’ ^32 equations (A29) to 

(A34) yields the following two simultaneous equations in terms of W]_]^, 
^21^ and ^31* 


1 + 




A2 ^(a3 - A2) 


2Al 


(w2i + JW51) 


A3/A2 


(A 35 ) 


and 


1 ± 




W21 + JW51 + — 


2Ai 


A2^(A3 - A2] 


1 ± 1 1 + 


^3^(^3 " ^2) / A2A3' 


2Ai 


exp \vjj_ 


Ai 


^21 " <3^31 2 


2Ai 


(A3 - A2) 


A2/A3 


(A36) 


Inasmuch as ^21^ and wjq represent the first switching point, 

these two equations express the locus of first switching points . 

Equations (A33) and (A3lf) can be manipulated to yield the wqW2 
and the W2W5 projections of the zero trajectory as follows: 

(,,2 . ^ 

A2^ 


exp 


^A2^A5 


Ai 


W12 + 


^2^3 1 

'^12; 


(A3T) 
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and 


/ 2 2\ " ^ 2 ) . p 1^2^^3 

(a^2 . ^^2) - ^ W22 = A32 exp\-^ ^^2 


A 2 M5" ^ 

^2 "^ 12 ; 


Al 


(A 58 ) 
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TABLE I.- TRANSFER-FUNCTION COEFFICIENTS 

b, radians/sec 6.65 

c, (radians /sec) 2 ^l 8 

Ki, deg/sec/volt l4.0 

K2, volts/g unit 0.517 

Kj, volts/g \mit/sec^ 0.47 

Ki^., deg/sec/deg 1,000 

g units/deg 905 

T, sec 0.1 
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Control-surface servo — ^ .< — Airframe — 



(a) Contactor roll control system. 


^ — Control— surface servo — ► ^ Airframe ^ 



(b) Optimum contactor roll control system. 

Figure 1.- Block diagrams of the contactor and the optimum contactor 
second-order roll control systems. 



Figure 2.- Block diagram of the optimum third-order contactor normaL 

acceleration control system. 
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Figure 3*- Typical response of the third-order contactor acceleration 
control system to a step command signal. 


ro 



Figure 4.- Block diagram of the limited-linear normal-acceleration control 

system. 


Q\l<L MI VOVM 


Error rate signal, e , deg/sec 


^1 £ 


— +1 

6 ^ 

^4 

a 

1 

J ^ 

.1 



S(TS+1) 










Error signal, deg 


Figure 5 *~ A phase-plane plot of the contactor roll control system in 
response to a 5 ° step input. The zero trajectory for this system is 
superimposed on the plane. 
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Time, t, sec 

Figure 6.- Transient response of the contactor roll control system to a 

5^ step input. 
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Figure 7*- Plot of the zero trajectories of the contactor roll control 

system on the phase plane. ^ 
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Figure 8.- Transient responses of the optimum contactor roll control system 

for step roll-angle Inputs. 
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Figure 9.- Initial switching curve in the plane of the optimum 

contactor normal— acceleration control system. 


38 


MCA TN 57^3 



-.03 -.02 -.01 0 .01 .02 .03 


Figure 10.- Initial switching curve in the WjW 2 plane for the optimum 
contactor normal-acceleration control system. 
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Figure 11.- Initial switching curve in the plane for the optimum 

contactor normal-acceleration control system. 
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W| W2 projection of locus 
of first switching points 


Locus of first 
switching points 


Figure 12.- Locus of first switching points in the upper half of the 
w phase space of the contactor noa^mal-acceleration control system. 



Figvire IJ.- Zero trajectory for the optinnim contactor normal-acceleration 

control system. 
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»2 

Figure l4.- A phase-plane plot illustrating the method of determining 
points on the first switching locus with an analog computer for the 
optimum contactor normal-acceleration control system. 
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Figure 15.- 


A typical trajectory in the plane for the optimum 

contactor normal-acceleration control system. 
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Figure l6.- A typical trajectory in the '^ 2^1 optimum 

contactor normal-acceleration control system. 
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Time^ t, sec 

Figure 17* - Normal-acceleration transient responses of the optimum con- 
tactor normal-acceleration control system to step inputs of n^. 
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Figure I8.- Normal-acceleration transient responses of the limited-linear 


control system. The gains 

:(t ) I dt 


mum value of 


K2 and were adjusted to yield a mini- 

for an input signal n^ of 
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